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ABSTRACT 

The  purpose  of  this  thesis  is  to  investigate  and  establish 
Sturm-Liouvil le  properties  for  special  eigenf unctions  which 
are  expressed  in  determinant  form.  In  particular,  a  special 
case  is  presented  where  the  elements  of  the  determinant  are 
Legendre  polynomials.  This  type  of  determinant  has  a 
probability  background  dealing  in  birth  and  death  processes. 
The  method  of  analysis  used  in  this  thesis  is  a  new  approach 
to  solving  this  specific  example.  This  investigation  involves 
systems  of  differential  equations  and  Prufer's  analysis  in  the 
phase  plane.  The  following  are  new  results  obtained  in 
addition  to  solving  the  special  case  mentioned  above. 

Special  determinants  of  hypergeometric  functions  also 
possess  Sturm-Liouvil le  properties.  As  a  special  case,  a 
different  proof  of  Turan's  Inequality  is  provided.  Finally, 
several  theorems  are  presented  for  Sturm-Liouvil le  systems  of 
differential  equations  with  polynomial  coefficients. 
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I .   INTRODUCTION 

A.   THE  PROBLEM 

The  main  purpose  of  this  thesis  is  to  investigate  and 
establish  the  essential  Sturm-Liouvil  le  properties  for  special 
eigenf unctions  which  are  expressed  in  determinant  form.  In 
particular,  a  special  case  is  presented  where  the  elements  of 
the  determinant  are  Legendre  polynomials  such  that 


^n,tix)- 


Pjx)    P„,iU) 


{n>k)       k=0,l,2,... 

.      (1.1) 
(-l<x<l)  n=l,2,3, . . . 


This  type  of  determinant  has  a  probability  background 
dealing  in  birth  and  death  processes  and  was  studied  by  Dr.  S. 
Karlin  and  Dr.  J.  L.  McGregor.  In  their  researches,  Karlin 
and  McGregor  encountered  determinants  of  classical  polynomials 
as  examples.  In  order  to  complete  their  results,  they  needed 
to  know  such  properties  as:  number  of  zeros,  interlacing 
zeros,  and  completeness.  Their  obvious  starting  point  was  to 
construct  an  ordinary  differential  equation  such  that  (1.1) 
was  the  solution.  However,  this  differential  equation  was  not 
in  a  form  that  suggested  how  to  answer  any  of  these  questions. 
[Ref.  1] 

The  method  of  analysis  used  in  this  thesis  is  a  new 
approach  to  solving  this  specific  example.  The  investigation 
that  follows  provides  a  constructive  answer,  not  only  to  this 

1 


example,  but,  to  a  whole  class  of  such  problems.  The 
following  are  new  results  obtained  in  addition  to  solving 
Karlin's  problem. 

Special  determinants  of  other  hypergeometric  functions 
also  possess  the  essential  Sturm-Liouvil le  properties. 

As  a  special  case  of  (1.1),  a  different  proof  of  Turan's 
Inequality  is  provided  in  Chapter  III.  Turan's  Inequality  is 
of  the  form 


P„,,  (x)  P„,.(x) 


<  0 


-l<x<l 
n=0 , 1/ 2, . . . 


(1.2) 


where   the   elements   of   the   determinant   are   Legendre 
polynomials. 

Finally,  several  theorems  are  presented  for 
Sturm-Liouvil  le  systems  of  differential  equations  with 
polynomial  coefficients.  In  general ,  boundary  value  problems 
of  systems  are  not  well-posed.  Until  now,  very  little  has 
been  established  concerning  these  systems  with  polynomial 
coefficients . 


B.   STRATEGY  FOR  SOLVING  THE  PROBLEM 

The  basic  strategy  developed  here  is  to  construct  a  second 
order  differential  equation  such  that  (1.1)  is  a  solution. 
Next,  this  differential  equation  is  reduced  to  a  special 
system  of  first  order  differential  equations  such  that  the 
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coefficients  are  polynomials.  Then,  the  solution  trajectories 
are  analyzed  in  the  corresponding  phase  plane  by  a  method 
known  as  the  "Priifer  analysis".  From  this  analysis,  several 
theorems  of  Sturm-Liouvi 1 le  type  are  developed  for  general 
systems  of  differential  equations  with  polynomial 
coefficients.  Finally,  completeness  and  orthogonality  are 
established  with  respect  to  a  positive  definite  matrix  weight 
function. 


II.   CONSTRUCTION  OP  A  STURM-LIOUVILLE  SYSTEM  OF 
DIFFERENTIAL  EQUATIONS 


A.   DEVELOPMENT  OF  A  NONSTANDARD  STURM-LIOUVILLE  DIFFERENTIAL 
EQUATION  WITH  SPECIAL  DETERMINANT  EIGENFUNCTIONS 

The  objective  of  this  section  is  to  construct  a  second 
order  linear  differential  equation,  with  appropriate  boundary 
conditions,  such  that  its  solution  is  (1.1).  For  the  general 
case,  with  arbitrary  n  and  k,  this  method  is  somewhat  lengthy 
and  requires  extensive  manipulation  of  equations.  However, 
considering  a  special  case  of  k  =  1,  this  method  is  not  as 
complicated  yet  still  illustrates  all  of  the  important  steps. 
For  completeness,  the  outline  of  the  solution  for  the  general 
case  the  same  procedure  and  is  provided  in  the  Appendix. 

In  the  general  case,  the  eigenf  unction  $jj]j(x)  has  the 
following  form 


*..ic(^) 


Ptix)   Pv,,(x) 


Pjx)    P„,,(x) 


n*l 


{n>k)       k=0,l,2,... 
(-l<x<l)  12=1,2,3,  .  .  . 


(2.1) 


Setting  k   =    1    results    in 


X         %(3^-l) 


♦«.i(^)  = 


Pnix)  P^I(X) 


The  strategy  is  to  use  recursion  formulas  and  to  develop 
three  equations  for  §jj|(x),  f'^^Cx),  and  i'^jCx)  as  functions 
of  Pj(x)  and  P'jj(x).  Then,  from  these  equations,  a  system  of 
equations  in  the  matrix  form  AX  =  0  is  constructed  where 
X  #  0.  Finally,  the  determinant  of  A  is  calculated  and  set 
equal  to  zero.  This  results  in  a  differential  equation  having 
(2.1)  as  a  solution. 

First,  the  equations  for  $j,j(x),  $'j,j(x),  and  $'jjj(x)  are 
developed.   Substituting  the  recursion  formula 


P,^,ix)=xPjx)-^.^^Fi{x) 


into  (2.1)  for  Pl,,(x)  and  P.,,(x)  results  in 


k+1 


•n,irU)  = 


PJX) 


D+1 


Pj,  ix)  xPj,  ix)  -  i^;^  Pi  (X) 


lE+r 


xP^(x)-il^Pi(x) 


Since  k  =  1,  P)j(x)  and  F' ^(x)  are  replaced  by  x  and  1, 
respectively.  In  addition,  by  computing  the  determinant  and 
collecting  like  terms,  we  obtain  the  following  equation  for 
$n,i(x). 


*«,l-  ^^'f^Pn(x)-     ""^llf^  P'n(x) 


n+l 


(2.2) 


Next,  we  differentiate  both  sides  of  (2.2).   Thence, 


<iU)  =  ^^'^^  pL(x)  -xP„U)  -il^PiU) 


(2.3) 


Then,  we  substitute  the  Legendre  differential  equation 


d 


il-jc')pi,{x) 


■-nin+l)P„{x) 


into  (2.3)  and  collect  like  terms.   Then,  (2.3)  becomes 


♦;,i(x)=(i3-i)xp,(x)  +  sb=iisii^p'^{x)     . 


(2.4) 


Repeating  the  same  procedure  on  (2.4)  results  in  the 
following  relation  for  ^"^^{k). 


<iU) 


PAx)    +  (n-l)xPi(x) 


(2.5) 


Now,  a  system  of  equations  in  matrix  form  is  constructed 
from  equations  (2.2),  (2.4),  and  (2.5).  The  system  takes  the 
form  AX  =  0 ,  i.e.. 


AX= 


1-x^ 


{n-l)x 
(n-l) (1-x^) 


2{n+l) 
-x(l-x^) 

J2+1 


-1  1 


PJx) 


pLix) 


=  0.  (2.6) 


Since  the  trivial  solution  can  not  be  a  solution  of  (2.6),  A 
is  not  invertible.  Therefore,  the  determinant  of  A  must  equal 
zero.  Hence,  by  setting  the  determinant  of  A  equal  to  zero  we 
obtain  a  second  order  differential  equation  where  i^^  ^(x)  is  a 


solution. 


\A\   =  {l-x^){l+2x^)^'iiix)    -   6Jc(l-x2)*i,iU)  + 


in-l) 


(3n+6)x^+i3 


-2kx 


(X)  =  0 


(2.7) 


By  setting  X   =  n-1,  then  (2.7)  takes  the  form 

(l-x»)  (1+3x2)  <i(x)  _  6x(l-x*)«i,i(x)  + 
X3(X+3)x2+X-lJ*„^i(x)  =0 


(2.8) 


Note  that,  for  each  Legendre  polynomial  solution  Pn(x)  to 
the  Legendre  differential  equation,  there  is  a  corresponding 
logarithmic  solution,  Qj|(x).  This  logarithmic  solution  obeys 
the  same  recursion  formulas  as  P^(x)  .  Hence,  by  the  same 
procedure,  the  second  solution  to  (2.8)  is 


Therefore,  the  general  solution  is 


Ci 


X      Vfe(3x2-l) 


Pjx)      P„,iU) 


+C, 


X      %(3x2-l) 


Onix)        J?^iU) 


(2.9) 


Finally,  translated  boundary  conditions  on  (2.9)  are 
imposed  so  the  second  logarithmic  solution  is  eliminated. 

Note  that  (2.8)  is  not  in  the  standard  Sturm-Liouvil  le 
form.  However,  {^^(x)  still  possesses  all  of  the  essential 
properties . 


B.   REDUCTION   TO   A   SYSTEM   OF   FIRST   ORDER   DIFFERENTIAL 
EQUATIONS 

The  objective  in  this  section  is  to  reduce  (2.8)  to  a 

system  of  differential  equations  of  the  special  form 


-1<JC<1 


(2.10) 


where  a,  IB,  and  y  ^^^  polynomials  in  x,  i  =  {^^(x),  and  T  is 
an  arbitrary  function  in  x.  There  are  an  infinite  number  of 
ways  to  transform  (2.8)  into  a  system  of  differential 
equations.   However,  by  requiring  polynomial  coefficients  we 
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remove  the  singularities  at  ±  (-3)"^.  Hence,  all  solutions  are 
well  behaved  on  the  interval  -l<x<l. 

The  strategy  is  to  convert  (2.10)  into  a  second  order 
differential  equation  and  then  equate  the  polynomial 
coefficients  of  i* ,  f'  and  §. 

The  system  (2.10)  has  the  following  equivalent  form 

(l-JC*)(3x*+l)»''-[(o+Y+2x)  (3jc2+1)+6x(1-x*)  1  ♦'  - 


*  =  0 


First,  note  that  the  coefficients  of  $'  are  the  same. 
Next,  equating  the  coefficients  of  $'  and  simplifying  results 
in 


Y  =  -2x-o 


(2.11) 


Equating  the  coefficients  of  §  and  substituting  (2.11)  for 
Y  leads  to 


3(X+3)x2+A-l 


-a'(3x2+i)  +a6x-jgiMil|i!? 
^    a(-2x-tt)  Ox'-H) 

1-JC* 


(2.12) 


Dividing  both  sides  by  3x'+l  and  simplifying  yields 


X[3(X-*-3)x^-^A-l]    _  /  ^  aSx 

(2.13) 
^(3x^-H)^  tt(-2JC-a) 
1-x^  l-^a 


Next,  we  consider  components  of  (2.13)  that  have  3x*+l  in 
the  denominator.   The  idea  is  to  select  a   so  that 

X[3(X+3)x^-*-k-l]-a6x   ■  0  modOx^+l) 

This  condition  on  o.  is  essential  in  removing  the  singularities 
at  ±(-3)"'.   For  now,  set  a  and  fl  so  that 

a  =px     ,    B=q 

where  p  and  q  are  constants.  Substituting  px,  p,  and  q  into 
(2.12)  for  Q,  a'  and  B,  respectively,  and  simplifying  results 
in 

(-3X2-9A.)x*+(2X2+10A.)jf2+(X2-X)- 

( -3p2-9p-9Jb)  Jf'+ ( -p2+2p-6Jb)  x^+ ( -p-jb) 

Next,  we  equate  the  coefficients.   Thus,  p  =  2X  and  q  = 
-A.(A.+1)   which  implies   that  a   =   2Xx  and  fl  =   -A.(A,+1). 
Substituting  back  into  (2.11),  it  follows  that  y  =  (A.+l)(-2x). 
Therefore,  the  system  (2.10)  takes  the  form 
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(1-x^)*'  =      2Xx^    +     (l+3x^)¥       -l<x<l 
(l-x^)f'  =  -X(X+1)*  +  (A+1)(-2a-)T 

A  simple  substitution  of  (A+l)f  =  f  implies  the  following 
equivalent  form 

(l-x^)*'  -  2kx^   +  (X+1)  (1+3x2) T      -i<x<l  . 

(2.14) 

(l-jc2)T'  =  -X*  +    (X+l)(-2x)¥ 

Hence,  (2.14)  is  a  system  of  first  order  differential 
equations  that  is  equivalent  to  (2.8). 

C.   PROPER  ANALYSIS 

The  Priifer  analysis  is  used  to  study  the  oscillatory 
behavior  of  trajectory  solutions  in  the  corresponding  phase 
plane  [Ref.  2:p.  312],  It  is  a  very  powerful  method  that  is 
extremely  useful .  The  objective  of  this  section  is  to  use 
this  method  to  establish  the  theorems  listed  below  for  general 
systems  of  differential  equations,  and  then  apply  these 
theorems  to  (2.14).  These  theorems  for  systems  are  adapted 
from  the  analogous  theorems  for  Sturm-Liouvil le  differential 
equations  [Ref.  2:pp.  312-353]. 
Theorem  1.   (STURM  SEPARATION  THEOREM) 

If  ($  f)^  and  (f  f)2  are  linearly  independent  solutions  of 
the  system 
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■--^it')  -  (I  ?)  (*)  -^<-<^ 


and 

4j8(-C>>-[Y-«]^>0 

where  a,   Q,   y ,    and  Q  are  polynomials  in  x,  then  i^  must  vanish 
at  one  point  between  any  two  successive  zeros  of  $2  •   ^^^  other 
words,  the  zeros  of  f^  ^^^^  *2  ii^terlace. 
Theorem  2.   (STURM  COMPARISON  THEOREM) 

Let  ($  T)  and  ((  j\)   be  nontrivial  solutions  of  the  systems 


-^>($;)-(i?)(?)--(-^)(3-(if)(^)-^^'<^ 


where  ol  ,    0,  y*  Q-  ot,  15",  Y'  si^i'^  Q  ^i^®  polynomials  in  x,  and 

4jg(-C>)-[Y-«]2>0 

4(jg-B)  (-(C>-^))-[Y-Y-(«-«)]*>0 


Then,  \    vanishes  at  least  once  between  any  two  zeros  of  f, 
unless  (^  T))  is  a  constant  multiple  of  (§  7). 
Theorem  3.   (OSCILLATION  THEOREM) 

Consider  the  associated  Priifer  differential  equation 
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(l-x')^=A.Scos*e+  [  (X+1)  y-Xa]  sinScose-  (A.+1)  ffsin'O 


obtained  from  the  system  of  differential  equations 


B(-^)($:)  -  {i%  it:ii?)  ({) 


(2.15) 


-l<x<l 


where  a,  fl,  Y'  ^^'^   Q  ^^^   polynomials  in  x.   If 

4JS(-C>)-[Y-a]»>0 

and  ©(-1)  =  X,  such  that  -%/2  <  x  <  %/2  for  each  X,  then  the 
solution  6  is  a  continuous  and  strictly  decreasing  function  of 
X   for  fixed  x  on  -l<x<l.   Moreover 


limeU).-  ^dlime(x) 


X-00 


for  -l<x<l. 

Theorem  4 .    Consider  the  system  of  differential  equations 

(2.15)  which  satisfies  the  Oscillation  Theorem.   If 

e(-l)-T,   -^<T<-^    for  all-^<X<» 
0(l)=o-i2ic,  -?<<'<^    I2=0,±l,±2,  .  .  . 


then  (2.15)  has  an  infinite  sequence  of  real  eigenvalues 
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. .  .<A..2<X.i<Ao<Xi<X2< 


with 


1^  ^«=-     and     ^^"^  ^»' 

J3-»oa  J3-»-« 


The  eigenf unction  $jj  belonging  to  the  eigenvalue  X^  has  exactly 
n  zeros  in  the  interval  -l<x<l,  and  is  uniquely  determined  up 
to  a  constant  factor. 
1.   Phase  Plane 

The  objective  of  this  section  is  to  convert  to  polar 
coordinates  the  following  general  system  of  differential 
equations 

(1-^).'  =    .*      .    c^    _^^^^^  (,,,, 

where  a,  fl,  Y'  ^^^   Q  ^^^    polynomials  in  x.   By  introducing 
polar  coordinates 

*  =r  ix)  cosO  ix) ,      T  =r  (x)  sinO  (x) 

where  6  is  positive  in  the  counterclockwise  direction,  the  new 
dependent  variables  become  r  and  0. 

Expressing  r  and  0  explicitly  as  functions  of  S  and  7 
leads  to 
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J 


ii) 


e-arctanf-Jl       and  (2.17) 


r2=»2+T2  .  (2.18) 

Then,  by  differentiating  both  sides  of  (2.17)  and 
(2.18),  an  equivalent  system  of  differential  equations  for 
(2.16)  is  presented  as  functions  of  r  and  0.  Differentiating 
both  sides  of  (2.17)  yields 

Since 


then  (2.19)  becomes 


dB  ,  W'^'^y  (2.20) 


By  substituting   from  the  general   system  (2.16), 
equation  (2.20)  takes  the  form 

CtX  ^3  +^^2 
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Substituting  in  polar  coordinates,  the  differential 
equation  for  8  becomes 

(l-xa)^=&os2e+(y-a)sindcos0-(?8in2O        ^^'^^^ 


Finally,  we  differentiate  both  sides  of  (2.18)  and 
follow  similar  steps.  The  differential  equation  for  r  then 
becomes 

(l-x*)^=r  [0008^6+  (CM-jg)  sinecose+ysin^e]    .    (2.22) 


Therefore,  (2.21)  and  (2.22)  represent  an  equivalent 
system  to  ( 2 . 16) . 

The  advantage  of  this  system  over  (2.16)  is  that 

(2.21)  is  a  first  order  differential  equation  of  only  0.  This 
allows  the  study  of  the  oscillatory  behavior  of  the  solution 
trajectories  independent  of  r.  To  find  r  for  any  solution  0, 
one  solves  (2.22)  directly. 

We  also  note  that  r  is  always  positive  for  nontrivial 
solutions.   This  can  be  seen  by  examining  the  solution  of 

(2.22)  [Ref.  2:p.  312]. 

Next,  we  recall  the  following  result  from  elementary 
algebra.   A  function 


F{x»y)    =  a  x^  *  b  xy  -*■  c  y* 


is  negative  definite  if  and  only  if 
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a<0  and    4ac  -  b^>0 

and  positive  definite  if  and  only  if 

a>0  and    Aac  -   jb*>0 

Next,  we  apply  this  result  to  (2.21).   If  we  can  show 
that 

4Jg(-C»-[Y-«]^>0 

then  equation  (2.21)  is  either  negative  definite  or  positive 
definite,  depending  on  the  sign  of  B.  If  equation  (2.21)  is 
positive  definite,  then  the  change  in  the  angle  6  is  always 
greater  than  zero.  This  means  that  the  angle  6  is  always 
increasing.  Since  our  polar  coordinate  system  is  set  up  so 
that  the  positive  direction  of  rotation  is  counterclockwise, 
then  the  solution  trajectory  of  6  always  rotates 
counterclockwise.  If  equation  (2.21)  is  negative  definite, 
then  the  solution  trajectory  of  B  always  rotates  clockwise. 
In  addition,  since  the  change  in  the  angle  ©  is  either 
positive  definite  or  negative  definite,  then  the  direction  of 
the  solution  trajectory  6  can  never  change.  Therefore,  we  can 
conclude  that  once  a  solution  trajectory  crosses  an  axis  it 
can  not  cross  back. 

Also,  since  f  =  rcos6,  then  the  zeros  of  §   occur  when 
6  =  x/2+n'K   where  n  =  0,  ±1,  ±2,  .... 
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Consider  the  system  of  differential  equations  (2.14). 
The  associated  Priifer  equation  is  negative  definite  for  X>0 , 
and  positive  definite  for  X<0.  However ,  for  (2.14),  it  is 
easy  to  determine  the  direction  of  rotation  across  the  S-axis 
and  f-axis  without  the  Priifer  equation.  First,  we  examine  the 
case  when  X>0  by  directly  substituting  into  (2.14). 

By  setting  $  =  0  we  can  determine  the  direction  of 
rotation  along  the  ?-axis.  Also,  by  setting  ¥  =  0  we  can 
determine  the  direction  of  rotation  along  the  f-axis.  Hence, 
if  we  substitute  these  values  back  into  (2.14),  then  the 
corresponding  phase  plane  for  A,>0  takes  the  form  in  Figure  1. 


¥ 
A 


M 


^  • 


Figure  1.    Phase  Plane  Analysis  Along  f-axis  and  f-axis 
for  A.>0 
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Clearly,  the  solution  trajectories  can  cross  the  f- 
axis  and  f-axis  only  in  a  clockwise  direction.  In  the  case 
when  A<0,  the  trajectories  rotate  in  the  counterclockwise 
direction. 

2.   Sturm  Separation  Theorem 

Suppose  (f  f)^  and  (f  f)2  are  two  linearly  independent 
solutions  of  the  system 

(1-;^)  (*'A    '   n   ?1  fS)   -1<X<1         (2.23) 


($;)  ■  (5  ?)  (?) 


and 

4jg(-C?)-[Y-«]*>0 

where  a,  6,  y,    and  Q  are  polynomials  in  x. 

There  are  two  cases  to  consider,  one  being  when  fl<0 
and  the  other  when  OO  .   First,  we  consider  the  case  when  Q<0 . 

We  examine  the  area  "A"  of  the  triangle  created  by 
joining  the  points  ($  f)j,  ($  f)2,  and  the  origin  in  the  ($,?) 
plane,  Figure  2. 

By  inspection,  the  area  "A"  is  equal  to  zero  only  when 
(S  f)^  and  (f  T)2  are  on  the  same  straight  line  through  the 
origin.  This  can  occur  with  ($  ¥)j  and  (f  T)2  both  in  the  same 
quadrant,  or  in  opposite  quadrants.  However,  this  can  never 
occur  since  ($  T)j  and  ($  f)2  are  linearly  independent,  so  "A" 
has  positive  area. 
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Figure  2.   Area  Principle  Analysis 

Recall  that  two  nontrivial  solutions  of  (2.23)  are 
linearly  dependent  if  and  only  if  they  are  located  on  the  same 
straight  line  through  the  origin  in  the  ($,f)  plane.  Hence, 
the  angle,  u,  between  the  linearly  independent  solutions, 
($  f)^  and  (§  T)2/  is  always  0<«<Jt.  In  other  words,  (f  T)j  and 
(f  f)2  are  always  in  the  same  rotating  half  plane  and  never 
cross  each  other. 

Analytically,  we  can  show  that  the  area  "A"  is  never 
zero  by  the  following  method.  Using  the  area  principle,  we 
calculate  the  area  "A"  using  the  following  formula 
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A=Vi 


1  *.  ▼, 

■*•     ^1     *  1 

♦i  ''i 

1  *,  T, 

=  % 

1                 X 

♦a   ^2 

10      0 

A                 m 

=  V4(«iT,-^,TJ 


(2.24) 


The  change  in  the  area  is 


dA 
dx 


=  y2[*lT2  +  •^▼i  -  *iTi  -  *2¥l] 


Substituting  in  f'^,   $'2'   '' i'  ^^^     ''2  ^^^"^     (2.16)   and 
simplifying  gives 


dA 
dx 


=Vi(an)  (♦lYj  -  ♦aYi) 


Substituting  into  this  equation  from  (2.24)  results  in 


-^=(«n)A 

dx 


Finally,  by  solving  for  "A",  the  following  formula  is  derived 
to  determine  the  area  "A"  at  any  point  -l<x<l; 


A'k  exp/f'(a+Y)cftl 


where  k  =  A(-l)  is  the  initial  area. 

Since  a  and  y   are  polynomials  and  -l<x<l,  the  area  is 
finite.   Moreover,  A(-l)  can  not  equal  zero  since  the  two 
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solutions  are  linearly  independent.   Hence,  the  area  "A"  can 
never  be  zero. 

Next,  we  recall  that  B<0  and 

Therefore,  d6/dx  is  negative  definite  which  implies  that  all 
solutions  rotate  in  a  clockwise  direction.  Also,  once  a 
solution  crosses  the  ¥-axis,  it  can  never  cross  back. 
Finally,  since  the  area  can  not  equal  zero,  the  angle  tt 
between  (f  T)^  and  (§  f)2  is  always  satisfies  0<o<x. 
Therefore,  as  the  solutions  rotate  clockwise  in  the  ($,f) 
plane,  (f  f)^  and  ($  f)2  always  cross  the  f-axis  in  the  same 
order.  Therefore,  the  zeros  of  f^  and  $2  interlace.  The 
analysis  for  the  case  when  fl>0  proceeds  in  exactly  the  same 
way  except  that  the  solutions  rotate  counterclockwise,  q.e.d. 
3.   Sturm  Comparison  Theorem 

Let  (i  f)  and  ({  t|)  be  nontrivial  solutions  of  the 
systems 

respectively,  where  a,  13,  y,  Q,  a,  "E,   Y/  and  Q  are  polynomials 
in  X,  and 
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4S(-C>)-[Y-«]*>0 
4(j8-S)  (-(C>-I?))-[Y-Y-(«-«)]^>0 

Again,  there  are  two  main  cases  to  consider.  One  case  is  when 
fi-F  >  0  and  the  other  when  B-U  <  0.  The  second  case  is  the 
same  as  the  first  except  that  the  solutions  rotate  in  the 
counterclockwise  direction.  Hence,  it  suffices  to  consider 
only  the  case  when  B-T3"  >  0  and  that  all  solutions  rotate 
clockwise.  The  objective  is  to  analyze  the  two  solutions  in 
the  ($,f)  plane  through  the  use  of  the  Priifer  equation. 

The  associated  Priifer  equations  for  (2.25)  are  listed 
below,  respectively. 

(l-x2)_l=jgcos2ei+(Y-a)sineiCosei  -  (Psin^e^    (2.26) 


(l-x^)^    -  Scos2e2+(Y-o)sine2COse2  -  ^in^e^   (2.27) 


By  subtracting  (2.27)  from  (2.26),  the  resulting 
equation  indicates  the  rate  of  change  of  the  angle  formed 
between  (S  f)  and  ((  r\) .  In  addition,  by  superimposing  one 
solution  on  top  of  the  other,  with  6j  =  %2  (Figure  3),  then  the 
equation  becomes 
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-  iO-'D)  sin^i 


The  resulting  equation  is  a  Priifer  equation.  This 
equation  facilitates  the  analysis  of  the  solution  vectors 
(f  T)  and  (^  r\)    by  examining  the  behavior  of  the  two  angles, 

Vi  and  O2  * 

Since  B-B"  >  0,  and 

4  (B-S)  ( -  iO-O) )  -  [y-Y-  («-a)  ]  2>0 

then,  d/dx(6p62)  >  0.  This  means  the  rate  of  change  of  the 
angle  between  (S  ¥)  and  (^  t|)  is  always  increasing,  and  ((  f)) 
"runs  away"  from  ($  f ) ,  Figure  4. 

Note  that  (f  ¥)  and  (^  i\)  are  solutions  of  two 
different  systems  of  differential  equations.  The  two  solution 
may  be  linearly  independent  or  dependent.  Hence,  one  solution 
could  "lap"  the  other  solution  one  or  more  times. 
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Figure  3.   Superimposition  of  Solutions 


>  ♦ 


Figure  4.  One  Solution  "Runs  Away"  from  the  Other 
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Recall  that  the  zeros  of  S  and  (  occur  when  (i  T)  and 
(^  t|)  respectively  cross  the  f-axis.  Therefore,  as  the 
solutions  continue  to  rotate,  between  any  two  zeros  of  i  there 
is  at  least  one  zero  of  i.  Note  that,  if  ($  f)  =  c(^  i\) ,  for 
all  -l<x<l,  then  the  zeros  occur  simultaneously.  Hence,  ^ 
vanishes  at  least  once  between  any  two  zeros  of  S,  unless 
((  1))  is  a  constant  multiple  of  (S  T) .  q.e.d. 
4.   Oscillation  Theorem 

We  consider  the  system  of  differential  equations  of 
the  form 


(2.28) 


<--">(;:) -(ti,'i:i!?)  (I)  — 

where   a,    B,    Y'    ^^^  Q  ^^^  polynomials   in   x,    and 

4j8(-C>)-[Y-a]*>0 


The  system  (2.28)  is  basically  the  same  as  (2.23) 
except  now  we  have  introduced  an  eigenvalue,  X.  The  objective 
now  is  to  vary  X  thereby  to  study  the  zeros  of  f .  The  Priifer 
equation    for    (2.28)    is 

(l-jfa)       1  =XjScos2e,  +  [  (X+l)  yXa]  sine.cose, 

dJf  (2.29) 

'  {X*l)  OBin'd^ 
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As  before,  the  zeros  occur  only  on  the  f-axis  which  is 
the  same  as  saying  the  zeros  occur  at  0  =  x/2+nx  for  n  =  0, 
11,  12,....  Next,  we  impose  endpoint  conditions  on  6  such 
that 


8(-l)-T,   -^<t<-5.   for   aIJ-«><X<«> 

^    ^  (2.30) 

6(1)  =o-i2«,  -:^<o<-|    i3=0,±l,±2,  . . . 

4U  4b 


so  that  the  angle  x  is  as  close  to  x/2  as  possible.  Note  that 
the  endpoint  conditions  are  computed  from  tanO  =  f/$. 
Now,  we  consider  a  second  system 


where  u  is  an  eignevalue.  The  Priifer  equation  for  this  system 
is 

(l-x*)  --7^=ji&7os*ea+  [  (ii+l)  Y-»io3  sine-cose, 

-(ji+DCteln^ea 

Suppose  that  x  is  fixed  and  the  differential  equation  (2.31) 
satisfies  the  following  endpoint  condition 

6(-l)-t,  -^<f<4  ^°^     -••<|i<w 
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Following  the  same  basic  proof  as  in  the  Sturm  Comparison 
Theorem,  we  subtract  (2.29)  from  (2.31),  and  equate  6j  to  O^ 
at  Oi(-l)  =  02(-l)  =  X.   Hence, 


(2.32) 


(l-x^)^(e,-e,)- 
(H-X)  [Scos^Q^*  (y-«)  sineiCOS0i-OBin^J 

Since 

then  the  discriminant  of  (2.32)  is 

(M-X)M4ig(-(?)-(Y-o)2]>0 


Therefore,  d/dx(^2~^0  ^^  either  positive  definite  or 
negative  definite.  Accordingly,  there  are  several  cases  to 
consider.  However,  by  examining  one  case  carefully,  all  other 
cases  follow  directly,  with  only  the  direction  of  rotation 
varying.   Thus,  it  suffices  to  show  only  one  case. 

Let  X  be  fixed.  Suppose  that  0<u<A  which  implies  that 
\x-k   <   0.   Now,  consider  the  case  when  fl<0. 

Since  fl<0  and  n-X<0,  then  d/dx(62-Oi)  >  0.   Hence,  by 
considering  the  conditions  that  0<vi  and  0<X,    it  is  implied 
that  all  solutions  rotate  in  the  clockwise  direction,  and 
(*  f)  "runs  away"  from  (^  r\)    as  X   goes  to  positive  infinity. 
Figure  5. 

28 


Figure  5.   (*  f)  "Runs  Away"  from  (^  r\) ,    X>0 

If  X<0 ,    then  ]x-X   >   0  which  implies  that 
d/dx(62~®i)  <  0-   Therefore,  when  X<0 ,    all  solutions  rotate  in 
the  counterclockwise  direction  and  (§  f)  "runs  away"  from  (^ 
r\)    as  X   goes  to  negative  infinity.  Figure  6. 

In  summary,  as  X  goes  from  -«  to  «,  the  angle  6^  goes 
from  •  to  -to  in  a  strictly  decreasing  manner.  Hence,  for 
fixed  X,  0  is  a  strictly  decreasing  function  of  the  variable 
X.   In  addition. 


lim  e»-« 

X-»oo 


and 


lim  6»<» 


q.e.d. 
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>  • 


Figure  6.   («  f)  "Runs  Away"  from  (4  r\) ,    X<0 

5.   Sequence  of  Eigenvalues  and  Interlacing  Zeros 

Consider  the  system  of  differential  equations  (2.15) 
which  satisfies  the  Oscillation  Theorem.  Now  we  continue  with 
our  analysis  and  maintain  the  conditions  of  the  previous 
proof.   Namely,  let  0<u<A  and 


-« 


0{-l)»t,  ^<t<-i,  foz  all   -«<>.<• 


BiD'o-im,    -?<<'<-f#  i3»0,±l,±2, . . . 


where  x  is  as  close  to  x/2  as  possible. 
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Now,  we  let  p  be  fixed  and  B<0.  Note  that  a  value  for 
X  which  satisfies  both  endpoint  conditions  is  an  eigenvalue 
for  the  system  of  differential  equations  (2.28).  Let  6j(-l) 
=  X.  Then,  as  X  increases,  the  solution  trajectory  of  8^ 
rotates  in  a  clockwise  direction.  Since  dO^/dx  <  0,  the 
trajectory  of  8j  eventually  crosses  6j(l)  =  o-nic  (n  =  0). 
Hence,  the  first  time  that  the  trajectory  of  8j  crosses  the 
line  ©j(l)  =  o-nx  is  when  n  =  0.  We  call  this  eigenvalue  X^, 
Therefore,  when  X  =  Xq  ,  the  solution  trajectory  of  8^  does  not 
cross  the  T-axis,  Hence,  there  are  no  interior  zeros  for  the 
corresponding  eigenf unction  $q  of  X^ ,    Figure  7. 


^(iH^-iht 


Figure  7.   Interior  Zeros  of  f 


0 
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Recall   that   for   each   X^,   the   corresponding 
eigenfunction  is  determined  by  the  formula 


♦^"rcos6 


Now,  as  X  continues  to  increase,  the  solution 
trajectory  of  6j  continues  to  rotate  in  the  clockwise 
direction.  Eventually,  the  trajectory  of  ©^  again  crosses  the 
line  6^(1)  =  o-njc  when  n  =  1,  or  0^(1)  =  o-x.  We  call  this 
eigenvalue  Xj .  Therefore,  when  X  -  Xj,  the  solution  trajectory 
of  ©2  crosses  the  f-axis  one  time.  Therefore,  the 
eigenf unction  S^  of  X^  has  one  interior  zero.  This  process 
continues  as  X  increases  to  «.  Similarly,  for  X<0,  the 
solution  trajectory  of  0^  rotates  counterclockwise  and  this 
process  continues  as  X  decreases  to  -•. 

In  summary,  as  X  goes  to  ±o»,  each  time  that  the 
solution  trajectory  of  0^  crosses  the  line  0^(1)  =  o+nic,  the 
corresponding  value  for  X  becomes  an  eigenvalue  of  the  system 
(2.28).  Hence,  there  is  an  infinite  sequence  of  Xj  for  n  =  0, 
±1,  ±2,...  for  which  the  second  endpoint  condition  is 
satisfied.  Therefore,  there  are  eigenvalues,  Xj,  such  that 
.  .  .  <X.^<Xq<X^<  .  .  .  ,  and  the  associated  eigenfunction  belonging 
to  Xj  has  exactly  n  interior  zeros  in  the  interval  -l<x<l. 
Moreover,  it  follows  that 

llmX.  «•     ,  limX_—«» 

"     and     ^       " 
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In  addition,  any  two  solutions  which  satisfy  the  same 
initial  conditions  are  linearly  dependent.  Therefore,  each 
eigenvalue  uniquely  determines  an  eigenf unction  up  to  a 
constant  factor.   This  completes  the  proof  of  Theorem  5. 

For  the  system  of  differential  equations  (2.14),  the 
discriminant  of  the  associated  Priifer  equation  is  positive. 
Therefore,  all  of  the  previous  theorems  apply. 

D.   ORTHOGONALITY  AND  COMPLETENESS 

The  objective  of  this  section  is  to  show  that  in  a  Hilbert 
space  the  eigenvectors  of  (2.14)  are  complete  and  orthogonal 
with  respect  to  a  weight  function.  [Ref.  3:pp.  344-353] 

1.   Orthogonality. 

The  system  of  differential  equations  (2.14)  is 
equivalent  to  the  following  form 


<--")(?;)  -  K"  '-1^)  (')  *   (o  '-1^)  (•)■  ''■'* 


For  convenience,  represent  (2.34)  with  the  notation 
(l-x«)Zl-XAZ;^+BZ;^  (2.35) 

where  Zj^  is  the  eigenvector   (f  T)   associated  with  the 
eigenvalue  X.      A  and  B  are  the  2x2  matrices  given  below 
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(?-  Kl^)     a^  B  -  (0  -:lf) 


In  addition,  we  let  L  stand  for  the  linear 
differential  operator,  such  that,  when  L  operates  on  Z  then  LZ 
is  defined  by 

LZ=  (l-x»)  Z'-XAZ-BZ=0 

The  general  strategy  is  to  find  the  adjoint  operator, 
Lt,  such  that  w'(LZ)  -  (L  w')Z  is  an  exact  derivative.  Then, 
we  integrate  both  sides  to  establish  orthogonality.  Finally, 
matrix  algebra  can  be  used  to  find  the  corresponding  weight 
function. 

Let  W  be  a  2  X  1  column  vector  such  that 

-  [  (l-x^)  W^] '-  [  (l-x2)  W"^  l-XW^^-WW* 

Then,  by  requiring  the  exact  derivative 

dx 

it  follows  that 
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Hence, 

For  general  ]x,    the  adjoint  differential  equation  is 
[(l-x*)HrT-  -^WA  -  wm  oz  (2.37) 


Next,  we  return  to  (2.36).   Substituting  w''(LZ)  and 
(LW)Z   into  the  equation  and  simplifying  yields 

ax 


which  simplifies  to 


ax 


Substituting  in  (2.35)  and  (2.37)  results  in 


(X-|i)  Rr»AZ=-^  [W^il-x^)  Z] 
dx 


Now,  we  integrate  both  sides.   This  leads  to  the  following 
orthogonality  relationship 

(2.38) 


iX-\i) p  W'JiZ'iW^{l-x')Z)  W   =0 
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If  the  integrand  is  well  defined,  then  the 
eigenvectors  W  are  orthogonal  to  those  of  Z,  with  respect  to 
the  matrix  weight  function  A.  Next,  we  determine  the  nature 
of  W. 

Expanding  (2.37)  leads  to 

il-x^)  t/=2xW-iiA'^-B^        .  (2.39) 

We  write  W  in  the  form  W  =  J£  where  £  is  a  2  x  1  column 
vector,  and  J  is  a  2  x  2  matrix  of  the  form 


'(-! !) 


Then,  (2.39)  becomes 

Then  we,  premultiply  both  sides  by  J  .   Note  that 
J'^  =  -J.   Thus, 

(1-x2)E'=2jcE-^hJA'JS+JB'"JE  (2.40) 


Since 


.'AV.(.;  I)   (?f  3^,V^)  (.?  I)   .   A 


and 
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JB^'J 


[O  0  J 


then  (2.40)  simplifies  to 

(l-x2)E'=|iAE+(B+4jd)E  (2.41) 

where  I  is  the  identity  matrix.   Next,  letting 

E-(l-jf*)-2Q      and 
E-4x(l-x*)  -3Q+  il-x^)  -2Q 

where  fi  is  a  2  x  1  column  vector  and  substituting  into  (2.41) 
gives 

(l-x^)Q'=»iAQ+BQ    .  (2.42) 

So,  Q  =  Z.  where  Z   is  the  eigenvector  corresponding  to  the 
eigenvalue  vl.      Back  substituting  results  in 

l^»  JE  -  J{l-x^)-^Q   -  J-(l-x*)-2Z^ 

If  A  #  y,  then  the  orthogonality  relationship  (2.38) 
becomes 

Since  A  #  u,  then  X-y  #  0  which  implies 
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When  A  =  u,  then 

'1  VT         J^A 


f     zj   ^^  ,  Z.  dx  #  0 


(2.43) 


since  the  integrand  is  always  positive.  In  addition,  since 
the  integrand  is  a  polynomial,  then  it  is  finite.  This 
implies  Z|  is  square-integrable  with  respect  to  the  matrix 
weight  function 


(l-x*)^ 


Initially,  the  (1-x*)"^  term  may  seem  to  imply  that  the 
integrand  is  not  well  defined  at  ±1.  However,  (1-x")  cancels 
in  (2.43),  and  all  the  remaining  entries  are  polynomials. 

Next,  we  determine  the  sequence  of  eigenvectors. 
Clearly,  X  =  1,  2,  3,...  all  yield  solutions.  This  is  based 
on  the  definition  of  X  such  that  X  =  n-1  where  n  =  2,  3, 
4,....   Now,  we  consider  X   =  -1,  -2,  -3....   Since 


for  the  Legendre  polynomials,  the  following  relationships 
hold. 
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"^l*!  "  ^-t(i*l)*ll   ■  ^- 


X-2 


(2.44) 


^X*2     "  ^-t(X*2)*lJ  ^-1-3 


(2.45) 


Substituting  (2.44)  and  (2.45)  into  the  original 
definition  of  i^  ^  yields 


♦n,l  = 


P^        P2 


Px*l    Px*2 


^1    P2 


P-X-2    P-X-3 


for  X   >    1.      Equivalently,  for  X  <  -l. 


« 


D.l 


Pi  P2 

PX'2       PX'3 


This  form  provides  the  solution  to  the  original 
differential  equation  (2.8).  Therefore,  there  exists  a 
sequence  of  eigenvectors 

•••'(*L (*).a  '(*L  '(*).  '(*)x  '(*l  '■■■'(*).  '■■■ 


where  ($  f)  has  polynomial  entries  of  degree  <  m+3. 
2.   Completeness 

Finally,  we  determine  that  the  sequence  of  continuous 
eigenvectors  (f  f)^  for  m  =  0,  ±1,  ±2,  ±3,...,  which  are 
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orthogonal  and  square-integrable  with  respect  to  a  matrix 
weight  function,  is  complete.  That  is,  we  show  that  every 
square-integrable  function  F  =  (1-x*)  (f  g)'  can  be  expanded 
into  an  infinite  series  such  that 


Recall  that  there  is  a  factor  (1-x*)  which  cancels  from  both 
sides  of  the  equation. 

Consider  the  column  vector  (p  q)'  where  p  and  q  are 
polynomials  of  degree  <  m+3.  There  exist  constants  c.^^  and  c^ 
such  that 


(5)  ■  [  -  (?)..  •  -  (?).  1 


is  a  column  vector  with  polynomial  entries  of  degree  <m+2. 
This  follows  by  equating  the  coefficients  of  X'  and  solving 
the  resulting  system  of  equations  for  the  coefficients  c.^  and 
C|.  Since  (i  f).,  and  (S  T)^  are  linearly  independent,  then  c.^ 
and  c^  exist.  By  mathematical  induction,  one  can  find 
constants  C|^  for  k  =  -m,  ...  ,  -1,  0,  1,  ...  ,  m  such  that 

(3-- (I).. -- (5).  -  (5)/M5);— .  (5). 
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Finally/  we  let  F  =  (1-x*)  (f  g)'  where  f  and  g  are 
any  continuous  functions  on  a  finite  closed  interval  -l<x<l, 
and  let  e^,  e^  ^  ^  •  Then,  by  the  Weierstrass  Approximation 
Theorem,  there  exist  polynomials  such  that 


(l-jc2)/-P  <  Ci  and 


il-x^)g-Q\   <  Cj 


for  all  X  on  [-1,1].   Hence, 


<-^'(l)-ra|-  (li:5J^:3hft)— <«^-«^'(i) 


which  implies 


(l-x')[g)    -  E   ^it  (!)     ^  »^  <«i'«a>  [\) 


Therefore,  F  is  approximated  uniformly,  within  an 
arbitrarily  small  distance  e  >  0,  by  a  linear  combination  of 
the  (f  f ) .  Hence,  the  sequence  {($f),}  for  m  =  ...  ,  -1,  0, 
1,  ...  is  a  uniformly  convergent  sequence,  which  implies  the 
sequence  is  mean  square  convergent  on  -l<x<l  with  respect  to 
an  integrable  positive  weight  function.  Recall  that  (1-x*)' 
cancels  in  the  integrand,  leaving  a  positive  definite  matrix. 
Therefore,  the  sequence  is  complete.   [Ref.  2:p.  353] 
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In  addition,  the  constants  C|^  are  determined  in  the 
following  way 

^1  f,T      J^A 


^*         -1  nT       CTK 


fzI-il^^Z^  dx 


42 


III.   RESULTS  AND  CONCLUSIONS 


A.   RESULTS 

We  have  shown  in  this  thesis  that  there  is  a  special 
function  which  takes  the  form 


*..i:U)  = 


P„{X)     P^iU) 


in>k)         k=0,l,2,... 
(-l<x<l)  n=l,2,3,  . . .   , 


(3.1) 


where  the  elements  of  the  determinant  are  Legendre 
polynomials,  and  possesses  the  essential  Sturm-Liouvil le 
properties.  In  addition,  several  interesting  results  are 
obtained  and  listed  below. 

1.   Hypergeometric  Functions 

The  Legendre  differential  equation  is  a  special  case 
of  the  hypergeometric  differential  equation.   This  suggests 
that  a  determinant  of  the  form  (3.1),  with  other 
hypergeometric  functions  as  elements,  might  possess  similar 
properties . 

Recall  that  the  recursion  formula 


P,,,(x)=xP,(x)-^^^Pi(A:) 


and  the  Legendre  differential  equation  are  used  extensively  in 
the  initial  development  of  (3.1).  In  fact,  it  can  be  shown 
that  since  the  hypergeometric  function  can  be  expressed  in  a 
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recursion  formula  of  the  form 

F(r+1  ;  x)  =  Air  :  x)F{z  -,  x)   -^  B{i  -,  x)F'{i   ;  x)    , 

then  the  determinant  (3.1),  with  other  hypergeometric 
functions  as  elements  and  appropriate  boundary  conditions, 
possesses  the  same  Sturm-Liouvil le  properties. 

2.  Turan's  Inequality 

The  proof  of  Turan's  inequality  is  a  straight  forward 
application  of  the  general  case  (3.1).  Setting  n  =  k+1  and 
substituting  back  into  (3.1)  results  in  a  determinant  that  is 
of  Turan's  form.  Recall,  from  the  previous  analysis  of  (3.1), 
that  there  exists  a  X  such  that  (3.1)  has  no  zeros  in  the 
interval  -l<x<l.   Hence,.  Turan's  inequality  is  established. 

3.  General  Systems  of  Differential  Equations 

In  general,  boundary  value  problems  of  systems  of 
differential  equations  are  not  well-posed.  Until  now,  very 
little  has  been  established  concerning  systems  with  polynomial 
coefficients.  In  Chapter  II,  several  theorems  are  developed 
and  proved  for  such  systems.  Hence,  if  the  conditions  of  the 
theorems  are  satisfied,  then  the  system  is  said  to  be  of 
Sturm-Liouvil le  form. 
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B.   CONCLUSION 

The  analysis  of  this  one  problem  presents  several  possible 
avenues  of  further  research.  Only  the  first  step  was  taken  in 
a  direction  which  could  result  in  establishing  several  new  and 
substantial  results.  For  example,  the  next  step  would  be  to 
establish  the  case  when  (3.1)  has  other  hypergeometric 
functions  as  elements.  Then,  one  could  continue  the  research 
on  Sturm-Liouvi 1 le  systems  of  differential  equations  by 
examining  the  possibility  of  expanding  the  necessary  and 
sufficient  conditions  for  such  problems.  Research  on  this 
subject  is  ongoing. 
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APPENDIX 
I .   PROOF  FOR  GENERAL  CASE  -  i^  ^ ( x ) 

A.   PURPOSE 

In  this  appendix  we  briefly  outline  the  procedure  for 
establishing  the  Sturm-Liouvil 1 e  properties  for  the  general 
case,  fjj|^(x).  The  basic  procedure  is  very  similar  to  that 
illustrated  in  Chapter  II.  Thus,  only  the  important  steps  are 
highl ighted. 


B.   DEVELOP  DIFFERENTIAL  EQUATION 

By  following  the  procedure  in  Chapter  II,  it  is  easy  to 
establish  the  following  relations  for  f^j^Cx),  f'jj]^(x),  and 


«  =  -(l-x2) 


n+1 


k+1 


I  -. 


(1) 


(72  +  1)  (ic+1) 


nr^k 


(2) 


*''=  (i2-ic) 


{r2-ic+l)  -,/ 


73+ 


+1)  pip  2x  p/p/   (Jc-iJ-H)  p  pi 

r~^"^*   (73  +  1)  (iC+1)  ^"^*     ic+1    ^"^* 


(3) 


Equations  (1),   (2),  and  (3)  lead  to  the  differential 
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equation 


(1-A«) 


P/+^|=|^(Pi)2  *''-(l-A«) 


l!l^^i*-(&<^*^' 


«^-t- 


(4) 


(fl-ic) 


(^-^.DP/^^P^,-  ^^-7^Vi)V^^  <^-)' 


♦  =0 


where  the  general  solution  is 


C, 


^it 

Pi:.l 

+ 

^2 

^. 

Pn^l 

•''jc   -Pjr*! 


Ofl   (?n*i 


Applying  appropriately  translated  boundary  conditions 
eliminates  all  logarithms,  Q^  ,  in  the  general  solution.  Now, 
setting  A,  =  n-k-1  and  letting 


1-jfi 


^I^TIP^^*^'^^*'   ^^ 


(5) 


ic+1 


TTi^*^^- 


(6) 


then  (4)  takes  the  form 


(1-A*)  •''-  (1-A«)  ^*'+  (X+1) 
1? 


X+ 


2Pi:(Pjt+-5^Pi)  1 


*=0 


(7) 
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C.   REDUCTION   TO   A   SYSTEM   OP   FIRST   ORDER   DIFFERENTIAL 
EQUATIONS 

Next,   we   reduce   (7)   to   a   system   of   first   order 

differential  equations  of  the  form 


(l-^)*'  =  a*  +  C>V     -l<x<l 


(8) 


where  a,  13,  Y'  ^^^  Q   ^^^   polynomials  in  x. 

First,  note  that  Q  and  Q'  have  no  zeros  in  common.  This 
implies  that  all  polynomial  coefficients  of  $,  $' ,  and  f'  are 
reduced  to  lowest  terms  and  only  simple  zeros  occur  in  the 
denominator.   The  proof  is  as  follows. 

If  Q' =0  ,  then  either  P/ =0  or 

-^Pi  +  P^  =  0    .  (9) 

If  Pjj'=0,  then  Q=0  only  when  P|j  =  0.  However,  it  is  easy  to  show 
that  Pjj  and  Pj^'  are  relatively  prime.   Hence,  when  Pj^'=0,  then 
Q  #  0. 

If  (9)  holds,  then  by  solving  for  P^  and  substituting  back 
into  Q  it  follows  that  Q=0  only  when  Pi5'=0.  But,  by  the  first 
case,  this  can  not  occur.  Therefore,  Q  and  Q'  have  no  zeros 
in  common. 

By  transforming  (8)  into  a  second  order  differential 
equation  and  equating  coefficients  of  f'  and  $'  it  is  easy  to 
establish  that 
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Y  =  -2x-a 


(10) 


Equating  the  coefficients  of  f  and  substituting  (10)  for  y, 
results  in  the  following  equation. 


/+  Q&   -  a{-2x-a)  _  acy 


X+1 


(11) 


-A(A+1)- 


^PtiPk'^^Pi) 


Now,  we  isolate  and  examine  the  components  of  the  equation 
that  have  Q  in  the  denominator.  The  idea  is  to  select  a  such 
that 


a(?'-(A+l) 


^Pt^PK^-S^Pi:) 


k+1 


=  0  modd?) 


Substituting  (6)  for  Q'  implies 


^*(l^^**^*)]"^^*^^  [2^*<^i^*l^^*^ 


ic+1 


=  0  modiQ) 


However , 


lc+1 


Pt*  Pt 


(12) 


and  Q  are  relatively  prime.   This  is  shown  by  the  following 
argument . 
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Set  (12)  equal  to  zero  and  then  solve  for  P^.  By 

substituting  back  into  Q  it  follows  that  (12)  and  Q  have  zeros 
in  conunon  when  P^  and  Pj^'  have  zeros  in  common.  However,  Pj^ 
and  P|j'  are  relatively  prime.  Hence,  (12)  and  Q  are 
relatively  prime. 

Therefore,  since  (12)  and  Q  are  relatively  prime,  then 
there  exists  a  polynomial  M  such  that 


Jc+1 


Pi- (X^l)  P^  -  Af|  ill^  (Pi)  «+pJ 


(13) 


By  examining  (13)  it  follows  that 


(X+l)Pj,+AfP/  =  0  mod  (Pi) 


Moreover,  since  Pj^  and  F^'    are  relatively  prime,  then  there 
exists  a  polynomial  N  such  that 


(X+1)    +  AfPj,  -  NPi 


(14) 


However,  by  the  Euclidean  Algorithm  for  Polynomials,  there 
exist  unique  polynomials  M^  and  N^  such  that 

M^P^   +  i^iPi  =  1    . 

From  (14),  this  implies  that 

itf«-(A.+l)Af,  and     ^=  (k-*^l)N^ 
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Hence,  there  are  unique  polynomials  M  and  N  that  determine 
a  such  that  all  the  zeros  of  Q  divide  out. 
Finally,  choose  B  to  satisfy  (11). 

D.   STURM-LIOUVILLE  PROPERTIES 

Although  it  is  somewhat  tedious,  it  can  be  shown  that 

Therefore,  all  of  the  theorems  previously  established 
apply.  In  addition,  the  properties  of  orthogonality  and 
completeness  follow  in  the  same  manner. 
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